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ABSTRACT 


This  investigation  is  a  study  of  the  problem  faced  by  a  circuit 
designer  who  in  the  analysis  and  synthesis  of  his  circuit  must  use 
parameters  which  vary  over  wide  ranges* 

A  technique  utilizing  a  digital  computer  resulting  in  design 
curves  is  developed.  A  second  technique  is  demonstrated  where  the 
same  design  curves  are  obtained  by  another  method  employing  an  analogue 
computer. 

Both  techniques  are  applied  to  a  neon  pi ot oconductor  circuit  and 
design  curves  are  obtained  with  each  method. 

The  writers  wish  to  express  their  appreciation  for  the  assistance 
and  encouragement  given  them  by  Professor  W.  Conley  Smith  of  the  U.  S. 
Naval  Postgraduate  School. 


ii 


TABLE  OF  CONTENTS 


Section  Title  Page 

1*  Introduction  1 

2.  Description  of  Digital  Computer  Design  Technique  3 

3.  Description  of  Neon-photoconductor  Circuit  10 

4.  Preliminary  Circuit  Analysis  22 

5.  Digital  Computer  Program  Technique  28 

6.  Interpretation  of  Digital  Computer  Data  42 

7.  Description  of  Electronic  Analogue  Computer  Technique  54 

8.  Conclusions  66 

9.  Bibliography  71 


lii 


U.'T  Of  I... LUSTRATIONS 


Figure  Page 

1.  Basic  Neon-photoeenductor  Circuit  10 

2.  Compensated  Neon- photoconductor  Circuit  11 

3.  Photoconductor  Resistance  Versus  Temperature  Curves  13 

4.  Photoconductor  Resistance  Versus  Neon  Current  Curves  15 

5.  Photoconductor  Resistance  Variation  Table  16 

6.  Comparison  of  Experimental  and  Approximate  Mathematical  18 

Curves  of  Photoconductor  Resistance  Versus  Neon  Current 

7.  Compensated  Neon-photoconductor  Circuit  22 

8.  Basic  Program  Flow  Chart  30 

9.  Basic  FORTRAN  Listing  "Program  Basic"  33 

10.  Sample  Digital  Computer  Printout  for  Basic  Program  34 

11.  Compensation  FORTRAN  Listing  "Program  Compa"  35 

12.  Main  Program  Flow  Chart  37 

13.  Hyperbolic  Approximation  FORTRAN  Listing  "Program  Neon"  38 

14.  Second  Order  Approximation  FORTRAN  Listing  "Program  Light"  39 

15.  Sample  Digital  Computer  Printout  for  Main  Programs  40 

16.  Sample  Digital  Computer  Printout  for  Programs  Basic  and  44 

Corapa 

17.  Design  Curve  for  Hyperbolic  Approximation  -  Compensated  45 

Curve 

18.  Design  Curve  for  Hyperbolic  Approximation  -  Uncompensated  46 

Curve 

19.  Basic  Neon-photoconductor  Circuit  43 

20.  Final  Design  Curve  for  Hyperbolic  Approximation  -  Digital  49 

Computer 

21.  Comparison  of  Resultant  Hyperbolic  Area  Plot  with  and  51 

without  Delta 

22.  Schematic  diagram  of  Analogue  Computer  Circuit  56 

iv 


23.  Table  of  Variable  Parameters  with  Associated  Potent io-  57 

meters. 

24.  Analogue  Computer  Output  Curve*  Z  s  .4  59 

25.  Analogue  Computer  Output  Curve ,  Z  s  1,0  60 

26.  Analogue  Computer  Output  Curve ,  Z  s  1.7  61 

27.  Analogue  Computer  Output  Curve ,  Z  s  2.5  62 

28.  Final  Design  Curve  from  Analogue  Computer  Solution  63 

29.  Comparison  of  Final  Design  Areas  from  Hyperbolic*  65 

Second  Order  and  Analogue  Computer  Solutions. 


v 


V^LL; 


vf  ^ 


11 

12 

13 

V 

Rl 

R2 

Rp  1 
Kp  2 
Z 
T 


tut as.1  series  current  in  ampet .„s 

neon  tube  branch  current  in  amperes 

shunt  branch  current  In  am pares 

maintaining  voltage  of  neon  tube  in  volts 

series  branch  resistance  in  ohms 

shunt  branch  resistance  in  ohms 

series  branch  photoconductor  resistance  in  ohms 

shunt  branch  photoconductor  resistance  In  ohms 

photoconductor  resistance  magnitude  multiplying  factor 

temperature  in  degrees  Centigrade 


1,  Introduet  ion 


System  reliability  is  the  keynote  in  this  dynamic  ag®  of  mechaniz¬ 
ation  where  unfailing  performance  is  expected  and  required.  However,, 
eomponent  parts  in  a  systems  configuration  are  often  subject  to  manufactur¬ 
ing  tolerances  which  may  vary  radically  from  the  norm,  yet  a  systems  de¬ 
signer  has  to  deal  with  these  varying  components  and  produce  a  dependable 
system.  Dependability  can  be  measured  as  a  statistical  quantity  or  as  an 
absolute  requirement.  The  latter  definition  of  dependability  will  be  used 
in  this  paper  as  a  yardstick  in  measuring  the  success  of  a  particular  de¬ 
sign.  The  problem  to  fee  undertaken  in  this  thesis  is  the  synthesis  of  a 
reliable  system  consisting  of  widely  varying  components.  Reliability  as 
used  here  is  dependent  on  the  components  themselves  not  failing. 

Circuit  synthesis  is  generally  interpreted  as  the  formulation  of  sig¬ 
nal  flow  charts  and  the  determination  of  component  values  for  the  realiza¬ 
tion  of  a  given  transfer  or  imaittanc©  function.  This  paper  will  be  con¬ 
cerned  with  only  the  latter  part  of  this  definition. 

The  specific  system  to  be  investigated  and  designed  consists  of  a 
neon  photoconductor  logic  circuit  which  is  a  basic  building  block  for  other 
logic  circuitry.  Justification  of  the  important  requirement  of  reliability 
in  this  circuit  is  indicated  by  its  use  as  a  basic  element  in  the  founda¬ 
tion  ©f  an  entire  system.  Analysis  and  synthesis  of  the  circuit  will  be 
restricted  to  a  steady  state  investigation. 

Using  the  specific  neon  photoconduetor  logic  circuit  as  an  example, 
a  general  technique  of  circuit  synthesis  will  be  proposed  in  which  the 
tools  of  circuit  analysis  and  inductive  reasoning  are  utilized  for  synthesiz¬ 
ing.  Analysis  can  exist  without  synthesis  but  the  converse  of  this  state¬ 
ment  is  not  true.  H.  S.  Scheffler  and  F.  R.  Terry  in  their  paper  entitled 
"Description  and  Comparison  of  Five  Computer  Methods  of  Circuit  Analysis" 
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describe  three  non- statist leal  approaches  to  circuit  analysis.  Some  of 
their  concepts  in  circuit  analysis  are  adopted  and  extended  into  the  field 
of  circuit  synthesis.  For  example  the  criteria,  point  of  failure,  is  a 
basic  extreme  condition  applicable  to  a  circuit  problem  whether  the  investi¬ 
gation  be  of  analysis  or  snythesis.  The  particular  logic  circuit  to  be 
investigated  will  be  synthesized  with  the  aid  of  both  a  digital  and  an 
electronic  analogue  computer.  Each  computer  has  advantages  over  the  other 
which  are  peculiar  to  the  methods  utilized  in  the  solution.  Successful 
circuit  synthesis  by  the  computers  produces  design  curves,  which  when 
properly  Interpreted  can  be  an  invaluable  asset  to  the  circuit  designer. 
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2®  Description  of  Digital  Computer  Design  Technique® 

Thif  section  explains,  t«  general  tenas,  a  technique  useful  in  cir¬ 
cuit  design  which  uses  to  full  advantage  the  high  speed  calculating  feature 
of  a  digital  computer® 

In  designing  a  circuit  to  meet  output  specifications  the  designer  is 
almost  invariably  faced  with  the  problem  of  using  components  in  the  cir¬ 
cuit  which  do  not  have  constant  values*  but  may  have  values  in  a  given 
range,  Th©  simplest  example  of  this  point  is  an  ordinary  fixed  resistance 
which  may  have  any  value  in  a  range  above  and  below  its  nominal  value® 

This  variation  is  due  to  manufacturing  differences®  A  range  of  variation 
may  also  be  due  to  environmental  changes  during  circuit  operation*  aging 
or  other  factors® 

Circuit  parameters  which  must  be  selected  from  within  a  known  range 
of  variation  will  be  referred  to  as  "known"  parameters®  Other  parameters* 
whose  value  can  be  chosen  by  the  designer  based  on  his  design  to  meet 
specifications  are  defined  as  "design"  parameters® 

The  final  circuit  or  circuits  decided  upon  by  the  designer  must  meet 
the  specifications  for  any  random  set  of  the  "known"  parameters  -  each 
parameter  in  th®  set  picked  arbitrarily  from  anywhere  in  its  range  of 
variation.  The  designer  must  choose  th®  design  parameters  for  his  circuit 
so  that  when  "acceptable"  values  of  the  design  parameters  &r®  used  in  the 
circuit,  the  circuit  performance  will  meet  specifications  for  any  random 
set  of  th©  known  parameters. 

In  simple  circuits  the  decision  as  to  acceptable  values  of  the  design 
parameters  may  be  easy.  Consider  a  simple  voltage  source  and  resistance 
circuit®  If  the  voltage  is  a  known  parameter  with  variation  from  5  to  10 
volts  and  the  specification  on  output  demands  a  current  of  fror.  3  to  4 
amperes,  then  the  selection  ©f  acceptable  values  of  resistance,  the  design 
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If  the  design  parameter,  R,  Tu-s  any  acceptable  value;  that  is  am- 
value  from  1.25  to  3.33  ohms,  the  circuit  will  meet  performance  specifica¬ 
tions  for  any  value  of  the  known  parameter,  voltage,  in  its  range  oc  5  to 

10  volts. 

The  example  shows  the  significance  of  the  terras  known  and  design 
parameters,  hut  is  much  too  simple  to  demonstrate  the  technique  developed. 
In  the  general  case,  there  are  many  known  parameters,  each  with  a  range  of 
variation.  There  is  usually  a  number  of  output  specifications  to  be  met 
and  the  designer  is  at  liberty  to  pick  his  design  parameters  in  sets  of 
acceptable  values. 

Stated  most  simply,  the  technique  begins  by  first  writing  equations 
which  express  the  output  as  a  function  of  the  design  and  known  parameters. 
By  inspecting  these  equations  and  employing  boundary  conditions,  the  de¬ 
signer  must  establish  a  range  of  investigation  for  each  of  the  design 
parameters.  This  range  is  similar  to  the  range  of  variation  of  the  known 
parameters,  but  differs  in  that  it  may  be  made  larger  or  smaller  during 
the  process  of  design  based  upon  the  designers  interpretation  of  the  in¬ 
itial  answers. 

Once  a  range  has  been  established  for  each  parameter,  the  equation 
expressing  the  output  must  be  solved  for  each  value  in  the  range  of  each 
parameter.  The  values  of  design  parameters  which  yield  acceptable  output 
regardless  of  the  values  of  the  known  parameters  used  are  recorded  and 
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plctte  1  *s  acceptable  t'onbinatt'Uja  '?  cN  ees i&n  parameter*..  k  circuit 
built  wfth  an  acceptable  c«rafe.ln^tl©n  of  these  design  parattiecar?"  will  pro- 
vide  output  which  meets  specifications  no  matter  what  values  of  known 
parameters  are  used. 


Since  it  is  Impossible  to  solve  the  circuit  equation  for  every  value 
©£  every  design  parameter  and  every  known  parameter,  the  ranges  of 
parameters  twist  be  divided  into  increments  and  the  equation  solved  for 
each  increment  of  each  parameter*  The  selection  of  these  increments  in¬ 
troduces  some  error*  As  in  other  approximation  procedures,  smaller  In¬ 
crements  introduce  less  error.  However,  smaller  increments  lead  to  more 
combinations  to  investigate  and  thus  many  more  calculations.  If  is 

the  number  ©f  increments  of  the  n  tb  parameter  to  be  investigated,  the 

* 

total  number  of  answers  to  b@  computed  and  tested  is 


It  is  seen  that  ®ny  reasonably  complicated  design  problem  with  a 
practical  number  of  increments  quickly  leads  to  an  extrema ly  large  nw.be r 
of  calculations ,  This  problem  suggests  use  of  a  digital  computer  with 
its  rapid  calculating  ability  as  an  invaluable  toole 

In  some  cases  it  may  not  be  necessary  to  investigate  all  increments 
of  a  known  parameter.  This  exception  was  demonstrated  in  the  trivial  ex¬ 
ample  stated  above  although  not  specifically  pointed  out.  The  range  of 
acceptable  resistance  was  determined  by  investigating  only  the  maximum 
and  minimum  values  of  the  known  parameter,  voltage.  Because  of  the  simple 
nature  of  the  circuit  it  is  obvious  that  a  value  of  the  design  parameter 
which  is  acceptable  at  both  the  maximum  and  minimum  values  of  the  known 
parameter  will  be  acceptable  for  any  value  between  and  only  the  end  points 
need  be  used.  The  characteristic  of  the  parameter  voltage  which  makes 
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is  always  of  the  same  sign. 


this  true  is  that  the  partial  derivative 


5 


Tits  ger^iat  then.  It  the  partial  derivative  of  the  output  taken  with 
respect  to  a  known  parameter  Is  always  of  the  same  sign*  for  all  values 
of  all  the  other  known  and  design  parameters,  only  the  maximum  and  mini- 
mum  values  of  that  known  parameter  need  be  investigated,  since  if  the 
design  parameters  are  selected  to  be  acceptable  with  both  these  points, 
they  will  be  acceptable  for  any  value  of  this  parameter  in  its  range. 

Since  parameters  of  this  sort  mean  a  great  reduction  in  the  number 
of  calculations,  the  equations  for  output  should  be  partially  different¬ 
iated  with  respect  to  each  known  parameter  in  turn  to  determine  if  deriva¬ 
tives  of  constant  sign  exist.  The  partial  derivative  must  be  of  the  same 
sign  for  al 1  values  of  the  other  known  parameters.  This  may  sometimes  be 
shown  by  manipulation  of  the  derivative  expression  to  show  it  is  always 
of  one  sign  due  to  its  form.  In  other  cases  a  worst  case  technique  will 
be  sufficient  to  prove  a  constant  sign  of  the  derivative.  In  this  pro¬ 
cedure,  if  the  derivative  is  assumed  to  be  always  positive  all  known  para¬ 
meters  appearing  in  the  expression  which  tend  to  increase  the  derivative 
are  set  at  their  minimum  and  all  those  which  tend  to  decrease  the  deriva¬ 
tive  are  set  at  their  maximum.  If  the  derivative  remains  positive  for  this 
worst  case,  it  can  be  assumed  always  positive.  This  procedure  has  the 
weakness  that  it  may  not  be  evident  whether  a  given  parameter  tends  to 
increase  cr  decrease  the  derivative.  Also,  in  general,  if  the  range  of 
parameter  variation  is  very  large,  the  derivatives  will  not  pass  this 
worst  case  test  and  the  parameter  must  be  divided  into  increments  for  in¬ 
vestigation. 

When  the  ranges  of  all  parameters  have  been  divided  into  increments 
to  be  tested  whether  numerous  increments  or  only  one  (end  points  type  of 
parameter)  the  computer  program  is  written.  This  program  solves  for  the 
output  using  each  increment  of  each  design  parameter  in  combination  with 
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each  lie  rement  ot  each  known  parameter.  The  output  resulting  from  each 
set  of  parameters  Is  tested  and  those  sets  which  result  in  output  which 
meets  specifications  are  stored  and  printed. 

The  output  of  the  computer  can  most  easily  be  interpreted  if  pre¬ 
sented  in  graphical  form.  If  there  is  just  one  design  parameter  this 
is  not  necessary.  A  most  useful  presentation  for  two  design  parameters 
is  a  plot  of  one  versus  the  other.  Starting  with  the  data  for  one  value 
of  a  selected  known  parameter;  for  each  value  of  one  design  parameter,  the 
maximum  and  minimum  acceptable  values  of  the  other  are  plotted.  Connect¬ 
ing  all  such  points  will  yield  an  area.  The  values  of  design  parameters 
corresponding  to  any  point  within  this  area  when  used  with  the  previously 
selected  known  parameter  will  result  in  a  circuit  which  will  meet  design 
specifications. 

This  procedure  is  repeated  until  area  plots  are  obtained  for  every 
increment  of  each  of  the  known  parameters.  In  each  case  one  design  para¬ 
meter  is  plotted  versus  the  other  utilizing  the  same  scales.  If  all  these 
plots  are  overlayed,  the  area  which  is  common  to  the  acceptable  area  of 
all  the  graphs  is  the  final  design  area  desired.  If  a  point  in  this  final 
common  area  is  selected  and  the  corresponding  values  of  design  parameters 
are  used  in  a  circuit,  the  circuit  will  meet  specifications  for  all  known 
parameters. 

The  requirement  for  obtaining  an  area  graph  for  each  increment  of 
each  known  parameter  demands  much  tedious  plotting.  This  is  essentially 
the  procedure;  however,  it  is  possible  to  write  the  computer  program  so 
that  the  data  will  be  partially  processed  by  the  computer  before  printing 
thus  eliminating  some  of  the  plotting  and  overlaying.  This  programming 
method  is  most  easily  described  by  use  of  a  specific  example  and  Is  de¬ 
monstrated  in  a  later  chapter. 
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It  .9  quit*  possible  that  when  the  acceptable  area  graphs  for  the 
increments  of,  known  parameters  are  over  1  ayed  there  will  be  no  acceptable 
area  which  is  common  to  all  the  graphs.  Thus  there  are  no  combinations 
of  design  parameters  which  will  guarantee  satisfactory  circuits  for  ran¬ 
dom  known  parameters.  This  condition  can  occur  because  the  specifications 
on  the  output  are  too  strict.  This  may  be  stated  another  way  by  saying 
that  the  range  of  one  or  more  of  the  known  parameters  is  too  large.  The 
difficulty  may  be  corrected  by  relaxing  the  specifications  on  the  output 
or  dividing  one  of  the  known  parameters  into  two  groups.  If  the  trouble¬ 
some  known  parameter  is  divided  into  two  groups,  each  of  which  has  half 
the  range  of  variation  of  the  parameter,  there  may  then  be  common  areas 
of  acceptable  values  for  these  half  ranges  and  all  other  known  parameters. 
In  this  case  then  there  are  in  effect  two  design  areas  rather  than  one  and 
the  divided  parameter  must  be  tested  and  separated  into  two  groups  before 
one  can  use  this  procedure  to  design  a  circuit. 

If  the  area  is  large,  it  may  be  possible  to  tighten  the  output  speci¬ 
fications  to  be  met  or  to  add  another  specification.  Doing  either  of 
these  will  tend  to  reduce  the  size  of  the  area.  By  inspecting  the  size 
of  the  area  resulting  from  the  Initial  output  specifications,  a  decision 
can  be  made  as  to  how  the  specifications  can  be  altered.  By  inspecting 
subsequent  results  and  making  further  adjustments  a  balance  between  specif¬ 
ications  and  available  combinations  of  design  parameters  can  be  achieved. 

The  technique  described  results  in  an  area  of  acceptable  design 
values.  If  the  designer  picks  a  set  of  design  values  from  in  the  final 
area,  he  is  assured  that  the  circuit  -silt  meet  output  specifications. 

The  question  still  remains,  "If  all  the  combinations  in  the  area  meet 
specifications  is  there  not  one  best  combination  or  some  part  of  the  area 
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which  yields  *£ttei  ijtcuits  chan  otner  parts?'*  The  techniques  des¬ 
cribed  up  to  the  point  will  not  answer  this  question. 

One  method  of  determining  the  best  part  of  the  final  common  area  is 
to  select  one  of  the  preliminary  area  plots.  This  plot  corresponds  to 
one  set  of  known  parameters.  Using  this  set  of  known  parameters  and  re* 
presentative  values  of  design  parameter  combinations  throughout  the  area, 
one  calculates  the  value  ot  one  of  the  output  quantities.  The  values  cat 
culated  are  plotted  at  the  corresponding  point  in  the  area.  Inspection  of 
the  plot  may  reveal  a  trend  in  the  value  of  the  output  parameter  and  in¬ 
dicate  a  portion  of  the  area  where  this  value  best  approximates  its  opti¬ 
mum  value.  This  process  may  be  repeated  with  other  sets  of  known  para¬ 
meters  and  the  corresponding  area  until  it  is  established  that  the  trend 
holds  in  enough  of  these  preliminary  areas  that  it  may  be  assumed  to  hold 
in  the  final  area  plot. 

In  summary  the  technique  which  has  been  described  in  this  section 
results  in  a  plot  of  one  design  parameter  versus  the  other.  An  area  of 
acceptable  combinations  of  the  two  design  parameters  can  be  drawn  on  this 
plot.  If  a  combination  is  selected  from  inside  this  area  and  a  circuit 
constructed  using  this  combination  and  any  combination  of  the  known  para* 
meters  used  in  the  circuit,  the  resulting  circuit  will  meet  the  desired 
output  spec i f ica t ions . 
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3.  Description  of  Heon-photeeonductat  Circuit. 

The  circuit  chosen  to  illustrate  the  technique  described  briefly  in 
the  preceding  chapter  is  a  neon  tube  *  photoconductor  circuit  which  is 
used  as  a  basic  building  block  in  the  design  of  computer  logic  circuits. 
The  basic  circuit  is  shown  in  figure  1* 


r 
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Figure  1. 

The  circuit  consists  of  a  neon  tube,  two  fixed  resistances*  R1  and 
R2  and  a  photoconductor,  Rpl,  supplied  from  a  300  volt  source. 

In  the  circuit,  the  output  quantity  of  interest  is  the  light  of  neon 
tube  two.  The  input  is  light  from  a  similar  external  neon  tube*  called 
number  one.  With  no  input,  there  is  no  output;  the  very  high  dark  resist¬ 
ance  of  the  photoconductor  which  is  in  the  mega  ohm  region  limits  current 
flow  and  keeps  the  voltage  drop  across  neon  number  two  low  enough  to  pre¬ 
vent  its  firing.  With  input  light  from  neon  number  one*  the  resistance 
of  the  photoconductor  drops  radically  to  its  light  resistance  value  of  a 


few  ktlohms  and  neon  number  two  fires,  producing  light  output « 

The  light  output  from  neon  number  two  is  used  as  input  to  sub- 
se  #jent  circuits  so  it  is  essential  that  if  be  as  nearly  constant  as 
possible  for  stable  operation*  This  constant  current  characteristic 
is  not  possible  with  the  circuit  as  shown  since  the  resistance  of  the 
photoconductor  is  a  function  of  temperature.  As  the  circuit  operates, 
the  ambient  temperature  rises*  causing  the  photoconductor  resistance  to 
rise  and  the  neon  current  to  decrease,  A  proposed  methcd  of  compensation 
intended  to  minimize  this  temperature  effect  is  to  insert  a  similar  photo 
conductor  in  the  shunt  branch  of  the  circuit  as  shown  in  figure  2* 


T 


?  a  *\r 


IV 


r 


i 


Figure  2. 

In  this  compensated  circuit  the  output  current  is  more  nearly  con* 
stant  despite  temperature  increase  because  of  the  feedback  nature  of  the 
output  light  shining  on  photoconductor  Rp2c  Now*  as  the  output  current 
tends  to  decrease  as  a  result  of  temperature  variation  changing  the  re* 
sistance  of  Rpl,  the  light  output  decreases,  increasing  the  resistance  of 
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Rp2  and  shunt  Ing  ot.  the  aval  able  current  through  the  neon  tending 

to  maintain  the  current  constant. 


The  problem  then  Is  to  design  a  circuit,  using  available  components, 
which  will  give  a  nearly  constant  output  light  or  current,  throughout 
the  operating  temperature  range  of  25  to  55  degrees  Centigrade.  The 
supply  voltage  is  assumed  to  be  a  constant  300  volts.  The  neon  to  be 
used  is  a  General  Electric  NE2  type  which  has  a  firing  voltage  of  135  volts 
and  a  maintaining  voltage  of  from  60  to  100  volts. 

The  characteristics  of  the  photoconductors  employed  in  the  circuit 
are  more  complex  since  the  value  of  photoconductor  resistance  is  a  func¬ 
tion  of  both  the  ambient  temperature  and  the  incident  light.  The  latter 
variation  is  more  easily  investigated  if  neon  current  rather  than  light 
is  considered  the  determining  factor.  Mathematical  relationships  governing 
photoconductor  resistance  variation  with  temperature  and  with  cucrent  have 
been  determined  from  experimental  data. 

The  temperature  dependence  of  resistance  can  be  obtained  by  consider¬ 
ing  figure  3  where  log^  of  resistance  in  kilohtifl  is  plotted  versus  tempera¬ 
ture  in  degrees  Centigrade  for  three  typical  photo  conductors.  The  curves 
are  all  for  a  constant  incident  light  corresponding  to  two  milliamperes  of 

neon  current.  Since  the  data  plots  as  a  straight  line  on  semi- log  paper 

r,  A  r~^1 

the  variation  is  of  the  form  K  p'  '  c-  .  If  curve  II  is  selected  to 
determine  the  constant. 


“f  A  /  #  ■'  * 


(1) 


(2) 


Dividing  (2)  by  (1) 


:K  a 


C  /  f  '/ 


Thus 


c  j 


12 


KiloMMS 
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The  coefficient  A  and  the  mathematical  relationship  expressing  re¬ 
sistance  variation  with  neon  current  can  be  obtained  by  considering 
figure  4.  In  this  eurve9  photoconductor  resistance  in  kilohms  is  plot¬ 
ted  versus  current  in  the  exciting  neon  expressed  in  milliampares ;  the 
curves  are  all  for  a  corstant  temperature  of  40°  C.  The  three  curves 
A,  8,  and  C  are  for  t^ree  representative  photoconductors ,  A  and  C  show 
the  extremes  of  variation  and  curve  3  is  considered  typical.  In  attempt¬ 
ing  to  express  the  photoconductor  resistance  as  a  function  of  light*  it 
was  noted  that  curves  A*  B  and  C  are  all  of  the  same  mathematical  form 
and  are  related  by  a  constant.  This  magnitude  multiplying  factor  will 
be  defined  as  Z.  This  relationship  is  best  illustrated  by  table  I,  figure 
5. 

The  current-resistance  characteristic  of  any  photoconductor  can  be 
obtained  by  deriving  a  mathematical  expression  which  will  fit  one  of  the 
curves  giving  the  proper  form  of  variation  and  then  multiplying  the  en¬ 
tire  expression  by  a  constant  to  fix  its  magnitude. 

Curve  B,  selected  as  typical  will  be  the  one  used  to  determine  the 
mathematical  form  of  variation.  If  this  expression  is  assumed  to  have  a 
magnitude  multiplying  factor  of  unity,  reference  to  figure  4  shows  that 
in  order  to  cover  the  range  from  curve  A  to  curve  C  the  magnitude  multiply¬ 
ing  factor  must  have  values  of  from  .4,  for  curve  A  to  2.5  for  curve  C  , 

The  first  approximation  to  curve  B  which  was  investigated  is  a 
hyperbolic  function  where  is  a  constant  to  be  determined. 

m  -  1  i 

^  }j  T 

Since  a  neon  current  of  2  ma.  will  be  used  as  the  mean  point  about  which 

to  design,  the  hyperbola  constant  was  determined  to  fit  curve  B  exactly 

3 

at  that  point.  From  tie  curve  at  I  ■  2,0  maM  Rp  ®  4.388  x  10  ohms. 
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Table  1 


Current 
in  nia. 

Rp 

Curve  A 
in  Kohm 

Rp 

Curve  B 

in  Kohm 

Rp 

Curve  C 
in  Kohm 

MA> 

Rp  (B) 

Rp  (C) 

Rp  (B) 

1.00 

4.48 

12.66 

30.8 

.354 

2.43 

1.25 

3.38 

9.07 

22.0 

.370 

2.43 

1.50 

2.78 

6.71 

17.0 

.414 

2.53 

1.75 

2.06 

5.32 

13.86 

.388 

2.60 

2.00 

1.72 

4.388 

11.22 

.394 

2.54 

2.25 

1.46 

3.79 

9.28 

.396 

2.44 

2.50 

1.26 

3.37 

7.88 

.374 

2.34 

2.75 

1.11 

3.04 

6.76 

.363 

2,22 

3.00 

1.00 

2.89 

5.82 

.346 

2.02 

f^i 

-  (?,I  .  i 

Figure  5 
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Incorporating  the  effect  of  temperature, 

£  7  74  K*  °  1x3  7 
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Since  the  value  of  R  as  4.388  x  10  ohms  from  curve  B  corresponds  to 


I  *  2ma.  and  T  ~  40°C  the  constant  8.776  can  be  determined 


,  3  /  j  x  /  ~ 


?.  77  L  At  - c  '  (  *  1 

fx  i  3 

?  77  c  A  2  *•  -X 

The  hyperbolic  expression  with  temperature  effect  included  approximat ing 


curve  B  is  then 


i ?, 


4  ?£  ,  -  -  j  *  1 

— _  c> 

X 

To  include  the  fact  that  the  magnitude  of  the  resistance  can  vary 
from  curve  A  values  to  curv®  C  values,  the  magnitude  multiplying  factor 


16 


Z  mast  be  Included.  The  complete  tesist&nce  expression  for  the  photo¬ 
conductor  Is  then 


In  the  circuit  being  analyzed,  Rpl  is  supplied  with  exciting  light 
from  an  external  neon,  neon  number  1,  which  is  assumed  to  provide  a  con¬ 
stant  light  intensity  corresponding  to  a  neon  current  of  2  ma.  There¬ 
fore:  y  -  ' 

j  l  A 


I 

* 


/ 


c  /  f  ? 


The  compensating  photo  conductor  Rp2  is  excited  by  light  fed  back 
from  the  output  neon,  neon  number  2.  This  light  corresponds  to  the  out* 
put  current  •  Therefore* 

y  2  '  c  /,j  y  / 

J  '?  _•  ±  £  2  ^ 

a/7  ^  - * - 

X  zz 

The  plot  of  the  hyperbolic  approximation  to  curve  B  shown  in  Figure 
6  indicates  that  the  hyperbolic  curve  fits  experimental  data  exactly 
only  at  I  «  2.0  ma.  Since  circuit  performance  with  currents  from  lc5 
ma.  to  2.5  ma.  is  of  interest,  the  use  of  the  hyperbolic  approximation 
introduces  some  error  especially  at  the  lowest  value  of  1.5  ma.  where 
the  error  is  12.8%. 

The  expression  for  photoconductor  resistance  is  used  in  subsequent 
calculations  and  the  simplicity  of  the  hyperbolic  expression  is  an  ad¬ 
vantage  when  used  in  more  complex  expressions.  However,  a  more  accurate 
curve  fit  over  the  entire  1.5  to  2.5  ma.  current  range  was  desired. 

The  second  approximation  to  curve  B  which  was  investigated  was  de¬ 
rived  by  assuming  the  curve  to  be  of  the  forms 

A  (?r  1 1  X  9  X  f  0  Zp  1-  p  -  ■ 

where  A,B,C,  and  D  are  constants  to  be  determined.  Appendix  I  shows  the 


1? 
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c&lcul&t  ioiois  Ics^iing  fco  th<s  final  £,:  cti 

,  V  ‘  .  ‘  ' 

from  which  f 

i  * 

L , 

j  ""  •  %  '  .  i  ~ 

Including  the  effect  of  temperature  as  in  the  hyperbolic  approxi¬ 
mation:  / 

/  .  (  5  ^  * 

l'  H  V  /  | 

^  :  1  9 

In  the  circuit  being  analyzed,  R1  is  supplied  by  light  from  neon 
which  has  a  constant  current  of  2.0  ma.  therefore 


resistance  have  such  large  ranges  of  variation,  an  output  current  of  from 
1.5  to  2,5  ma9  will  meet  specifications.  R1  and  R2  nrnst  be  chosen  such 
that  this  specification  is  wet  with  any  neon  tube,  and  any  photoconductor 
having  characteristics  within  the  ranges  discussed  above.  The  only  assump¬ 
tion  being  that  both  photoconductors  have  the  same  value  of  magnitude  multi* 
plying  factor,  Z.  A  further  specif Lea t ion  is  that  for  any  finished  circuit 
which  will  have  one  value  of  neon  voltage,  V,  and  one  value  of  Z,  the  vari¬ 
ation  in  current  due  to  temperature  change  from  25°C  to  55°C  shall  not  be 
greater  than  .15  ma0  The  effect  of  imposing  stronger  restrictions  on  the 
output  will  be  investigated  and  the  results  compared. 

The  upper  limit  on  the  values  of  R1  to  be  investigated  can  be  estab¬ 
lished  by  setting  the  current  II  *  1.5  and  the  voltage  V  *  60, 
then 

O,  '  •*.  .  , 

K  i  y  —■ 

It  .  X  « 

R1  .  cannot  be  as  conveniently  set  since  its  value  depends  upon  the  value 
min 

of  R2  which  has  not  yet  been  determined.  A  rough  initial  estimate  can  be 
obtained  by  assuming  no  current  flow  through  the  shunt  branch  and  ignoring 
the  contribution  of  Rpl  in  comparison  with  R1  then 


In  view  of  the  approximations  made  here,  a  more  pessimistic  value  of  20 

kilohms  should  probably  be  chosen.  Experience  has  shown  that  a  range  of 
3  3 

from  40  x  10  to  130  x  10  ohms  for  R1  will  adequately  cover  the  possible 
values  of  R1  to  have  the  circuit  meet  output  specif ications* 

The  maximum  value  of  R2  is  indeterminant  since  the  effect  of  increas¬ 
ing  R2  is  to  shunt  more  current  through  the  neon  tube.  It  can  be  seen  that 
if  II  is  2.5  ma,  or  less,  R2  can  be  increased  without  bound  and  the  current 
12  will  never  exceed  the  output  specification.  The  minimum  value  of  R2  is 
also  indeterminate ,  In  the  uncompensated  circuit  a  minimum  value  of  R2  can 


20 


be  determined  because  a  lower  limit  exists  below  which  the  neon  voltage 

will  not  he  above  the  135  volts  necessary  for  firing*  In  the  compensated 

circuit  however*  due  to  the  very  high  dark  resistance  of  the  compensation 

photoconductor,  the  neon  will  fire  for  any  value  of  R2.  The  Magnitude 

of  R1  ,  was  selected  as  an  initial  lower  limit  for  R2  also*  An  initial 
min 

upper  limit  of  700  K  ohms  was  selected*  Subsequent  experience  showed 

3  3 

that  a  range  of  R2  from  40  x  10  to  530  x  10  ohms  was  sufficient  to  yield 
the  desired  design  curves* 
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4 .  P  re  1  ii!T»iM  r  y  r:  i.  rc  u  1  Ami '  y  s  l  n . 

This  section  deals  with  the  application  of  the  design  technique  dave- 
loped  in  Section  2  to  the  me©n~photoeonduetor  circuit  which  has  been  pre~ 


viously  described* 

r 


Figure  7# 

The  output  quantity  of  interest  is  the  light  from  neon  number  2. 

As  previously  noted*  the  output  is  more  easily  expressed  as  the  current 
through  the  aeon,  an  equally  good  measure  of  performance.  The  technique 
first  requires  an  equation  expressing  the  output  as  a  function  of  the 
other  parameters. 
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If  the  hyperbolic  approximation  Is  used  for  photoeonduetoy  resistance,, 


JL  ^ 


h'i  *- 


>  '  y.  h  ?  /  / 
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This  equation  can  be  easily  manipulated  into  the  form  of  &  quadratic 

c'-  .[.  y'  i  /  C 


where 
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Then  by  the  quadratic  equation: 


«  o 


I 


“  1  j  _  }L  _ 

Inspection  of  the  coefficients  A»B,  and  C  shows  chat  since  A  is 
always  positive  and  C  is  always  negative  (the  maximum  value  of  V  is  100) 
the  quantity  resulting  fro®  evaluating  the  radical  will  be  greater  than 
B.  Since  current  12  cannot  physically  fee  in  the  direction  opposite  to 
that  assumed,  a  negative  answer  for  12  has  no  significance  and  the  aqua* 


tion  to  be  solved  with  the  computer  is 

t  z  +  fyy.AQ 

y  A 

If  the  second  order  approximation  for  photoeonductor  resistance  is 


used: 

Ij 
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This  equation  can  also  bs  manipulated  into  the  quadratic  form. 


•'  i  •:  s  :  '  ■ 

where  the  coefficients  are  somewhat  more  complex, 
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The  same  quadratic  formula  must  then  be  solved. 
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In  order  to  determine  whether  any  of  the  known  parameters  fall  ln~ 
to  the  category  where  only  the  maximum  and  minimum  values  need  to  be  in¬ 
vestigated,  partial  derivatives  of  12  with  respect  to  V,  Z,  and  T  must 
be  taken  and  cheeked  to  see  if  they  are  of  constant  sign* 
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For  the  hyperbolic  case: 
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Inspection  shows  chat  the  quantities  A*  1 and  — —  are  always 
positive*  C  and  -  jra  always  negative  and  the  sign  of  S  rmiy  be  plus 
or  minus  depending  upon  the  M.-gnii  ude  of  the  negative  term 
compared  with  the  por>trve  term 

t%  i!  \  x  i  «  •  ->•  ~  A  A  ;•  ✓  I 

The  large  number  of  terms  in  the  .  expression  make  it  unwieldly 

?  l 

and  no  manipulation  which  would  establish  the  sign  of  j  suggests 

4  7 

itself.  Therefore,  the  second  approach  must  be  used.  The  slope  will  be 
assumed  negative  and  the  "worst  case"  values  of  PI,  P.2,  T,  Z,  and  V  which 
will  tend  to  make,  the  slope  positive  will  be  inserted.  The  expression  with 
these  values  inserted  is: 
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The  expression  evaluated  is  a  +  .0374.  A  positive  answer  shows  that 
the  initial,  assumption  of  a  negative  slope  ir.ay  be  incorrect.  The  informa¬ 
tion  gained  from  this  test  is  not  conclusive,  since  it  is  improbable  that 
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the  worst  values  of  all  the  pu  rasters  will  occur  simultaneously. 

However*  since  a  positive  slope  was  calculated,  the  pessimistic  alter¬ 
native  must  be  chose*,  and  the  temperature  must  be  investigated  in  in¬ 
crements.  A  similar  worst  case  analysis  on  L  leads  to  the  same 

.>  X" 

conclusion. 

The  variation  of  X2  with  V  will  be  investigated  next. 
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tive  as  long  as 


are  always  greater  than  zero 
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Sjlii*.  e  A  ant  .•  i  f  hi <?  last  Inequality  holds  and  thus  the  in¬ 
itial  assumption  of  a  negative  -  is  true.  Since  '  is  always 

negative,  only  the  maximum  and  minimum  values  of  V  are  investigated* 

Wien  a  similar  analysis  of  these  partial  derivatives  is  made  using 
the  more  accurate  second  order  approximation,  it  is  found  that  , 

and  ;  all  have  possible  changes  in  sign  and 

J  >  J  ' 

therefore  v;hen  this  approximation  is  used  all  three  parameters  must  be 
investigated  using  increments.  This  condition  is  then  one  disadvantage 
to  using  the  second  order  approximation* 

When  it  has  been  determined  that  a  parameter  must  be  investigated 
using  increments,  the  magnitude  of  the  increments  must  be  chosen.  As 
pointed  out  in  Section  2,  the  number  of  increments  selected  is  a  com¬ 
promise  between  accuracy  and  keeping  the  number  of  calculations  reason¬ 
able.  The  initial  choice  of  increments  was  as  follows: 

T  =  25,  35,  45,  55, 

Z  ~  .4,  .9,  1.4,  1.9,  2.5, 

V  =  60,  100,  for  hyperbolic  approximation  , 

V  =  60,  70,  80,  90,  100,  for  second  order  approximation, 

Rl  a  in  10  kilohm  increments, 

R2  =  in  20  kilohm  increments* 

Plots  made  using  these  increments  showed  that  parts  of  the  range  of 
Rl  caused  the  value  of  R2  to  increase  rapidly  with  possible  inaccuracy  in 
the  plot.  The  reason  for  this  rapid  increase  is  explained  in  Section  6. 
To  eliminate  this  possible  inaccuracy  the  increments  of  Rl  were  altered 
so  that  in  the  questionable  region  the  increments  were  taken  much  smaller 
(2  Kilohm  increments). 

The  problem  has  at  this  point  been  reduced  to  the  solution  of  a 
quadratic  equation  with  coefficients  made  up  of  the  five  parameters  Rl, 
R2,  T,  Z  and  V,  The  equation  is  to  be  solved  for  all  combinations  of 
each  of  the  increments  of  the  five  parameters. 
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5,  Digital  Computer  Program  Technique 

The  solution  of  a  quadratic  equation  whose  coefficients  are  func¬ 
tions  of  five  variables  and  where  the  number  of  solutions  is  the  product 
of  the  parameter  increments  is  ideally  suited  to  a  digital  computer,  A 
parameter  increment  Is  defined  as  the  number  of  values  that  a  variable 
will  be  assigned  spanning  its  entire  range  of  variation.  The  number  of 
combinations  of  parameters  to  be  examined  in  this  case  will  be  denoted  by 
N  where: 


N 


J 


U  £  ' 


Programming  of  a  digital  computer  in  machine  language  requires  an 
elaborate  and  complex  sequence  of  instructions  to  specify  every  step  of 
an  operation.  With  the  advent  of  the  FORTRAN  processor  in  late  1956, 
the  intricacies  of  programming  were  no  longer  a  stumbling  block.  Al¬ 
though  the  FORTRAN  system  was  developed  by  J.  W.  Backus  at  IBM,  it  has 
found  universal  acceptance  for  all  types  of  digital  computers.  The  FORTRAN 
processor  has  made  it  possible  for  scientists  and  engineers  to  apply  the 
language  of  ordinary  mathematics  to  their  individual  problems. 

Quoting  from  a  paper  prepared  by  Backus  and  his  associates,  "The 
purpose  of  FORTRAN  was  to  reduce  by  a  large  factor  the  task  of  preparing 
scientific  problems  for  IBM’s  next  large  computer,  the  704.  If  it  were 
possible  for  the  704  to  code  problems  for  itself  and  produce  as  good  pro¬ 
grams  as  human  coders  (but  without  the  errors),  it  was  clear  that  large 
benefits  could  be  achieved.  For  it  was  known  that  about  2/3  of  the  cost 
of  solving  most  scientific  and  engineering  problems  on  large  computers 
was  that  of  problem  preparation.  Furthermore,  more  than  90  per  cent  of 
the  elapsed  time  for  a  problem  was  usually  devoted  to  planning,  writing 
and  debugging  the  program.  In  many  cases  the  development  of  a  general 
plan  for  solving  a  problem  was  a  small  job  in  comparison  to  the  task  of 
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devising  and  *;  •  .0  machine  procedures  to  carry  out  che  plan.  The  goal 

of  the  .ORTFAN  project  was  to  enable  the  programmer  to  specify  a  num¬ 
erical  procedure  using  a  concise  language  like  that  of  mathematics  and 
obtain  automatically  from  this  specification  an  efficient  704  prog 
to  carry  out  the  procedure.  It  was  expected  that  such  a  system  would 
reduce  the  coding  and  debugging  task  to  less  than  one- fifth  of  the  job 
it  had  been." 

FORTRAN  (derived  from  formula  translation)  is  a  language  that  great- 
ly  resembles  the  language  of  ordinary  mathematics.  Symbols  such  as, 

*,  />  +>  in  FORTRAN  are  digested  as  exponentiation,  multiplication, 
division,  addition,  and  subtraction.  Mathematical  operations  such  as 
square  roots,  logarithms,  trigonometric  and  hyperbolic  functions  are 
available  in  the  computers  library  of  routines  and  used  merely  by  re¬ 
ferencing  the  appropriate  operation  by  abbreviated  name. 

The  Control  Data  Corporation  1604  digital  computer  was  married  to 
the  FORTRAN  processor  in  late  1961.  The  results  of  the  merger  produced 
a  high  speed  computer  that  could  be  commanded  with  relative  ec.se.  Suc¬ 
cessful  results  with  FORTRAN  programming  can  best  be  achieved  by  follow¬ 
ing  these  basic  rules  in  the  order  shown. 

1.  Block  diagram  -  (flow  chart) 

2.  Mathematical  and  logical  operations  to  be  performed 

3.  Transformation  to  FORTRAN  language. 

An  example  of  a  basic  flow  chart  for  the  solution  of  the  quadratic 
equation  using  the  hyperbola  approximation  is  illustrated  in  figure  8. 

The  idea  behind  this  program  is  to  see  the  effect  on  the  neon  current  as 
the  temperature  changes  for  any  one  set  of  known  and  design  parameters. 
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Fig.  8 
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The  flow  chart  indicates  the  sequence  of  operations  to  be  performed 
with  alternate  paths  lependent  on  logical  tests*  Each  block  in  the  flow 
chart  may  represent  a  few  calculations  or  logical  interrogations.  Es¬ 
sentially  the  flow  chart  in  figure  8  systematically  takes  every  possible 
combination  of  circuit  parameters  and  tests  the  output  current  to  see  if 
it  meets  specif ication3 .  If  a  failure  occurs  the  design  parameter  is 
altered  and  the  calculations  and  logical  tests  are  started  over  again* 

When  the  entire  range  of  temperature  variation  has  been  investigated  and 
yielded  an  acceptable  circuit,  the  design  and  known  parameter  set  which 
was  responsible  for  the  acceptable  circuit  for  the  entire  temperature  varia 
tion  is  recorded  on  magnetic  tape®  In  addition  to  the  parameter  set  the 
variation  in  neon  current  is  also  placed  on  the  magnetic  tape.  This  pro¬ 
cess  continues  with  a  new  parameter  combination  until  all  parameter  combin¬ 
ations  have  been  investigated. 

The  mathematical  calculations  to  be  performed  in  the  solution  of  a 
quadratic  equation  consist  of  evaluating  the  coefficients  which  are  vari¬ 
ables  dependent  upon  the  parameter  set  under  consideration  and  then  using 
the  quadratic  formula  to  evaluate  the  neon  current.  Since  in  the  present 
problem  the  argument  of  the  square  root  is  always  positive  and  greater 
than  the  coefficient  B,  only  the  positive  root  need  be  calculated.  Logi¬ 
cal  tests  to  be  performed  consist  of  determining  acceptable  circuit  designs 
by  employing  the  given  circuit  specifications  as  test  criteria.  Other  log! 
cal  tests  are  used  in  determining  maximum  -  minimum  values  and  program  loop 
completions. 

The  transformation  from  mathematical  to  FORTRAN  statements  is  accom¬ 
plished  next.  Input  to  the  1604  digital  computer  can  be  by  either  punch¬ 
ed  cards  or  magnetic  tape.  Utilization  of  punched  cards  is  superior  at 
this  time  as  cards  which  are  in  error  can  be  easily  replaced*  and  addi¬ 
tions  or  deletions  made  as  necessary.  The  complete  program  in  FORTRAN 
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appears  as  in  figure  9*  bach  line  in  the  program  represents  a  single 
punched  IBM  card*  If  this  same  program  had  been  done  in  machine  lan* 
guage,  these  thirty  odd  FORTRAN  cards  would  be  replaced  by  about  five 
hundred  machine  language  cards. 

If  the  FORTRAN  program  when  placed  in  memory  is  complete  and  logi¬ 
cal,  the  program  will  compile.  If  compilation  is  not  accomplished  an 
error  printout  will  result*  The  error  printout  categorizes  the  errors 
that  have  escaped  detection  by  the  programmer  and  lists  these  errors  by 
FORTRAN  statement  number.  When  the  error  is  recognized  alterations  to 
the  punched  cards  can  easily  be  accomplished  and  the  program  compilation 
retried  until  successful*  Writing  output  on  magnetic  tape  is  a  separate 
entity  from  the  computer  proper  and  occurs  simultaneously  with  computer 
calculations.  The  magnetic  tape  in  turn  is  read  out  on  an  off  line 
printer  so  as  to  conserve  computer  time.  An  example  of  the  digital  com¬ 
puter  printout  as  obtained  from  the  off  line  printer  appears  as  in  figure 
10. 

prior  to  proceeding  with  these  digital  computer  answers  a  mathemati¬ 
cal  check  is  required  to  see  whether  the  problem  is  being  solved  correct¬ 
ly.  A  desk  calculator  should  be  used  rather  than  a  slide  rule  to  obtain 
enough  significant  figures  for  comparison  purposes.  When  a  calculator 
check  has  been  made  for  one  set  of  parameters  and  dependability  of  the 
computer  program  proven,  then  the  accuracy  of  the  remaining  output  can  be 
relied  upon.  A  similar  FORTRAN  program  for  the  uncompensated  neon  photo¬ 
conductor  circuit  is  illustrated  in  figure  11.  This  program  is  needed 
for  a  comparison  of  its  results  with  those  of  the  compensated  circuit. 
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PROGRAM  BASIC 

DIMENSION  VI  3)  ,R  I  (  10)  ,'ft2  I  1  7)  ,T(4  )  ,  Z(  4  ) 

100  FORMAT ( 3F4.0 ) 

READ  1 00  < V( I ) . 1= 1. 3) 

200  FORMAT  (  10F7.0)  ,,, , 

READ  200  ( R  1  ( J ) •  J=1,  lOT 
300  FORMAT  (  11F7.0) 

READ  300 ( R2( K 1 »K= 1 » 17) 

400  FORMAT  (4F3.0) 

READ  400  (TIL),  L  =  1 . 4  )' 

500  FORMAT  ( 4F3. 1 ) 

READ  500  { Z ( M  > , M= 1,4) 

DO  16  1=1,3 
DO  15  M=  1 , 4 
DO  14  J  =  1,10 
DO  13  K=1,  17 
CM  I N 1  =  3. 

CMAX  1=  1.. 

DO  12  L  =  1,4 

A  =  C  R 1 ( J )*R2l K 1+2. 1 1 E  3*  Z ( M )*R2(K)*E XPF 1 -0 1 83*  T ( L) 1 )/ 1 .E6 
E  =  2. 1 1 E3*Z ( M) *EXPF( .0183* TIL ) 1 »( V(  I )  M2.*R  1 1 J)  1/1.E31+VI I  1 »R 1 1 J) 

F  =  8.904E3*Z(M)*Z(M)*EXPFl .0366*T(L ) )+R2(K)*(V( I  1-300.) 

8  =  ( E  +  F  )/  1.E3 

C  =  4.22*Z(M)»EXPFI .0183*TIL) 1  *  I  V  ( 11-300. ) 

CTWO  =  ( -B+SQRTF ( ( B*B)-(4.*A«C) > 1/(2. *A) 

IF  (CTWO-1.5) 13, 13,8 

8  IF  (CTWO-2.5 19, 13,13 

9  CMAX  1=MAX 1 F ( CMAX 1 ,CTWO ) 

CMIN1=MIN1F(CMIN1,CTW0) 

12  CONTINUE 
DELTA=CMAX1-CMIN1 

600  F0RMAT(F5.0,2F13.0,F8.  1,F13.5) 

WRITE  OUTPUT  TAPE  4 ,600 , VI  I ) ,R 1 ( J ) ,R2 I k ) , Z I M 1 , DELTA 

13  CONTINUE 

14  CONTINUE 

15  CONTINUE 

16  CONTINUE 
END  FILE  4 
PAUSE  1 
END 

END 

60.  80.100. 

40000.  50000.  60000.  70000.  80000.  90000.100000.110000.120000.130000. 
200000.220000.240000.260000.280000.300000.  320000 . 340000. 360000 . 3BU000 . 400000 . 
420000.440000.460000.480000.500000.520000. 

25.35.45.55. 

2.51.71.0  .4 
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PROGRAM  COMPA 

DIMLNSION  V(  3  )  ,R.M.lO)  ,R2(17),T(4),7(4) 

100  FORMAT ( 3F4 .0 )  ‘  ’ 

READ  IOO(V{  I  ),  I=j1, 3) 

200  FORMAT  (  I0F7.0  ) 

READ  200  (  R  1  {  J  )  t  J=,1 ,  10) 

300  FORMAT  l  1  IF  7.0)  '  ' 

READ  300  { R2  (  K  ) ,  K=.l  »,1 7  ) 

400  FORMAT  (4F3.0) 

READ  400  {  T(  L  )  ,  L=  I  ,,.4.) 

500  FORMAT  ( 4F3.  1  )  ’ 

READ  500  (Z(M),M=.1,4) 

DO  16  1=1,3 
DO  15  M=  1, 4 
DO  14  J  =  1,10 
DO  13  K=  1 ,  I  7 
CM  I N 1  =  3. 

CM AX  1=  1  . 

DO  12  L  =  1,4 

OCTWO={ {300. -V{ I )  )/(Rl{J)+2.  I  I E  3#  Z  {  M) *EXPF(.01H3*T(L) ) )  —  ( V ( l )/ 

1  R2 { K  )  )  )  *  I  .E3 
IF  { CTWO—  1.5)  13,  13,0 

8  IF  ( CTWO-2 .5)9, 13,13 

9  CMAX1=MAX1F(CMAX I, CTWO) 

CMIN1=MIN1F(CMIN I, CTWO) 

12  CONTINUE 

CELT  A=CM AX 1-CMIN 1 

600  FORMAT (F5.0, 2F  I  3 . 0 , F8 . 1 , F I  3 . 5 ) 

WRITE  OUTPUT  TAPE  4 , 600 , V ( I  )  , R  1  ( J ) , R 2 { K ) , Z ( M ) , OE L T A 

13  CONTINUE 

14  CONTINUE 

15  CONTINUE 

16  CONTINUE 
END  FILE  4 
PAUSE  1 
END 

END 

60.  80.100. 

40000.  50000.  60000.  70000.  80000.  90000.100000.110000.120000.130000. 
200000.220000.240000.260000.280000.300000. 320000.340000.360000.380000.400000. 
420000.440000.460000.480000.500000.520000. 

2.51.71.0  .4 
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The  rain  c^mputa:  program-  illustrates  one  of  the  prime  i  r  mintages 
of  the  digital  computer  in  evaluating  circuits  which  have  known  para¬ 
meters  of  the  type  requiring  incremental  analysis.  The  flow  chart  of 
the  main  program  called  neon  or  light  which  uses  the  hyperbola  or  second 
order  approximation  is  illustrated  in  figure  12.  Here  again  every  possi¬ 
ble  combination  of  circuit  parameters  are  systematically  tested  to  deter¬ 
mine  if  they  meet  specif icat ions 0  With  each  set  of  design  parameters 
the  known  parameters  requiring  increments  are  combined  in  every  possible 
manner.  If  every  combination  of  these  incremental  known  parameters  pass¬ 
es  all  logic  tests  with  one  set  of  design  parameters,  then  this  combina¬ 
tion  of  design  parameters  is  recorded  on  magnetic  tape  along  with  the 
maximum  and  minimum  neon  current  for  the  particular  voltage.  This  signi¬ 
fies  that  the  present  design  parameters  being  considered  have  produced  a 
circuit  which  has  successfully  passed  all  of  the  logic  tests  for  the  un¬ 
predictable  known  parameters  Z  and  T.  The  FORTRAN  programs  for  neon  and 
light  are  depicted  in  figures  13  and  14. 

Each  set  of  output  data  in  figure  15  is  a  successful  design  combina¬ 
tion.  In  other  words  for  a  neon  tube  maintaining  voltage  of  60  volts, 
a  series  resistance  of  90  [<  ohms,  and  a  shunt  resistance  of  200  K  ohms,, 
the  neon  current  will  only  vary  between  2.32944  and  2,01816  milliamperes 
for  the  hyperbola  approximation  and  between  2.32938  and  2,01802  millamper- 
es  for  the  second  order  approximation*  This  variation  in  neon  current  in¬ 
cludes  any  variation  in  the  characteristics  of  the  photoconductor  due  to 
temperature  or  light  intensity  within  the  assumed  possible  limits.  The 
last  two  columns  of  the  output  data,  maximum  and  minimum  neon  currents  for 
the  chosen  design  parameters,  have  no  immediate  use  in  the  synthesis  of 
the  circuit.  Since  these  currents  had  to  be  calculated  to  evaluate  the 
success  of  the  circuit,  it  was  an  easy  task  to  store  them  in  memory  and 
print  out  the  maximum  and  minimum  values  for  each  design  parameter  set* 
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Fig.  12 
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PROGRAM  NEON 

DIMENSION  V  ( 3  ) ,  R  1  (  2  1 ) »  R2 ( 34 ) ,  Z  (  5  ) 
FORMAT ( 3F4 . 0 ) 

READ  100(V( I ), 1=1,3) 

enD  M  AT  I  1  IC7  n  \ 


4 

5 


FORMAT  (  1  1F7.0) 

READ  200  (R  I  ( J  )  , J= 1 ,21 ) 

READ  200  ( R2 ( K  )  , K= 1 , 34 ) 

FORMAT  ( 5F3. 1 ) 

READ  500  (Z(M),M=1,5) 

DO  16  1=1,3 
DO  15  J=  1 , 2 1 
DO  14  K=  1,34 
CM  I N  1  =  3. 

CM AX  1  =  1  . 

DO  13  M=l,5 

A1  =  (Rl(J)*R2(K)+2.  1  1E3*Z(M)»R2(K)*EXPF( .4575) )/l .E6 
E1=2. 1 1E3*Z ( M)*EXPF( .4575 )*( V(t)+( 2.*R1 ( J) ) /l.E3)+V( I) «R1 ( J) 
F1=8.904E3*Z(M)*Z(M)*EXPF(  .9150)  +R2  (  K  )  *  (  V  (  I  )-300.  ) 

B 1  =  ( E 1 +F  1  )  / 1 . E3 

Cl=4.22*Z(M)*EXPF(.4575)*(V(l )-300. ) 

CT1=(-B1+SQRTF( (B1*B1)-(4.*A1*C1)))/(2.«A1) 

I  F ( CT  1-  1  .5  )  14,4,4 

IFtCTl  —  2.5)  5.5,14 

A2=(R1(J)»R2(K)+2.  1  1E3»Z(M)»R2(K)*EXPF( .6405) )/l.E6 
E2=2. 1 1E3*Z(M)*EXPF( .6405) » (V(I)+(2.*R1(J))/1.E3)+V(I)»R1(J) 
F2  =  8.904E3*Z(M)*Z(M)*EXPF(  1.281)+R2(K)*(V(I )-300. ) 

B2= ( E2+F2 ) / 1 . E3 
C2  =  4.22»Z(M)*EXPF( .6405 )*( V( I ) -  300 . ) 

CT2=(-B2  +  SQRTF( ( B2*B2 )- ( 4 . *A2*C2 ) ) ) / ( 2.  *A2 ) 

IF ( CT2-  I  .5 )  14,6,6 

6  I F ( CT  2-2 . 5 )  7,7,14 

7  A3= ( R 1 ( J  )  *R2 ( K ) +2.  1  1E3»Z(M)»R2(K)*EXPF( .8235) )/ 1 .E6 

E3  =  2 . 1 1  E3*Z ( M )  *E  XPF ( .8235)*(V(I)+(2J*R1(J))/1.E3)+V(I)»R1(J) 
F3=8.904E3*Z(M)*Z(M)*EXPF( 1 . 647 ) +R2 ( K ) * ( V ( I )-300. ) 

B3= ( E3  +  F  3  )  /  1  . E  3 

C3  =  4.22*Z(M)»EXPF( .82  35 ) *( VI  I ) -300 - ) 

CT3=(-B3+SQRTF( ( 83* B3 ) - ( 4 . * A3*C3 ) ) ) / ( 2 . * A3 ) 

I  F ( CT 3-  1  .5 )  14,8,8 
IF ( CT3-2.5 )  9,9, 14 

A4=( R1  ( J  )  *R2(K  )+2. 1  1E3*Z(M)»R2(K)*EXPF (  1 .0065 ) ) / 1 . E6 
E4  =  2  .  1  1E3*Z(M)*EXPF(  1.0065)*(V(l)-+(2.*R1(J))/1.E3)  +  V(I)*Rl(J) 
F4=8.904E3»Z(M)*Z(M)*EXPF(2.013)+R2(K)*(V(I ) -300 . ) 

B4= ( E4  +  F4  )  / 1 . 1 3 

C4  =  4.22*Z(M)*EXPF( 1 . 0065  ) * ( V ( I )-300.  ) 

CT4=(-B4+S0RTF( ( B4* B4 )- ( 4 . * A4»C4 ) ) ) / ( 2 . * A4 ) 

IFICT4-1 .5)  14,10, 10 

10  I F ( CT4- 2 . 5 )  11,11,14 

11  CMAA=MAX1F( CT 1 ,CT2»CT3,CT4 ) 

CM  I N=M I N 1 F ( C T 1 , C T2 , CT3 , C  T4 ) 

DELTA=CMAA-CMIN 
IFIDELTA-. 15)  12,12,14 

rMAVl  —  MAVICfTMAA  TMAVl 
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13 
600 

14 

15 

16 


I F ( DELTA- . 1 5 )  12,12,14 

CM AX  1=MAX 1 F { CMAA,CMAX 1  ) 
CMIN1=MIN1F(CMIN,CMIN1 ) 
CONTINUE 


CONTINUE 

FORMAT ( F6.0, 2F 1 3.0,  2F 1 3.  5 ) 

WRITE  OUTPUT  TAPE  4,600,  V  (  I')„R  1  (  J  ) ,  R2  (  K  )  ,CMAX  1 ,  CMI  N 1 
CONTINUE 


WRITE 
CONTINUE 
CONTINUE 
CONT  INUE 
END  FILE  4 
PAUSE  1 
END 
END 

60.  80.100. 

40000.  50000.  60000.  70000.  75000^.  80000.  82000.  84000.  86000.  88000. 
92000.  94000.  96000.  98000.  100000*.  1  05000.  1  1 0000.  1  1  5000. 1 20000 . 1  30000 . 
40000.  60000.  80000.  100000.  120000. 140000.  160000.  180000.200000.220000. 
260000.280000.300000.320000.340000.360000. 380000.400000.420000.440000. 
480000 .500000 .520000. 540000. 560000'.  580000.  600000.620000.640000.660000. 
700000. 

.4  .91.41.92.5 
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PROGRAM  UGH  I 


100 

200 

500 


4 

5 


0  I  MfcNS  I  ON  V(  5  )  ,R  1(  2  I  )  ,  R2  (  7,4  )  ,ZC  5  ) 

F  0  R  M  A  I  (5F4.0) 

Rt:  AU  100  (  V  (  I  )  •  1=1,5) 

FORMAT  (  1  IF/.O ) 

RFAO  200  ( R  I  ( J  )  , J=  1 ,2  I  ) 

REAU  20U  (R2CK)  ,K=1 ,3M 

FORMAT  (5F3.II 

RF AD  500  (Z(M),M=I ,5) 

DO  16  1=1,5 
DO  lb  J=  1 , 2  I 
DO  1  4  K=  1 , 34 
CM  IN  1  =  3. 

CM  AX  1=1. 

DO  13  M=l,5 

OA 1=2. 1 I  * Z  (  M ) *  E  XP F ( . 45  7  5  )  * (  1 . 42*R2 ( K) +  5 3 0. 2 « / ( M ) * FXPF ( .45  75 ) + .25 1 18 
1  * R 1 ( J ) )  +  1.42*K1(J)*R2(K)/l.t  4 

EI=2795.6*ZCM)»EXPF(. 4  575)*(R1CJ)/I.E 3  +  2. 1I«ZCM)*EXPF(. 45751-1300. 
1  —  V (  I  )  )/5.273)-RI(J)*R2(K)/1.E3 

Fl=2.  I  1*ZCM)«EXPF(  .  4  5  7  5  )  *  (  1  .  4  2  F  3  *  V  (  I  )  - R  2 ( K )  )  +  I  . 4  2  * ( V (  I  )  »  R  1  ( J  )  -  4  0  0 . 
1 *R2 ( K  )  +V (  I  ) » R 2 ( K  )  ) 

R1=E1+F  1 

C  1  =  7.  (  M)  *EXPF(  .4578)  *(  2  795.6*  (VC  I  1-300. 1-2.  1  1k3  + 

I  VC  I )  )  +  ( 300.- VC  I )  )  *  R  2 ( K  )  -  V ( I  )  *  R  I  (J) 

CT  1=  ( -H  1  +  SQR  TF  (  (  11  l*R  1  )-(  4.*  A  1»C  1  )  )  )/  (2.+A1  1 
I r ( C  T  1—  1 . 5  )  14,4,4 

IFCCT1-2.5)  5,5,14 

A2  =  2.  I 1»Z(M)*FXPF( .6405 (*(  1 . 42 *R2 ( K) +530. 2  * / ( K ) * EXPF l . 6405 ) + .25 1  1 8 
1 *R  1  ( J  )  )  +  1 .42*R  1  (  J  )  *R2<  K  )  /  1  .1'  4 

E 2  =  2  795. 6* Z ( M)*FXPF( .6405 ) *(K  I  ( J ) / 1.F3  +  2. 11 *Z CM) *EXPF (. 6405 )-( 300. 
1  — V ( 1  )  )/5.273)-Rl(J)*R2(K)/l. E3 

F2  =  2. 1 1 »Z ( M ) *EXPFC .6405  )  * (  1 .42E 3*V ( I  )-R2(K)  )  +  1 .42# ( VC l ) *R1 ( J )-300. 
1  * R 2 ( K ) +  V ( I ) «R2 ( K  )  ) 

R2=F2+F2 


10 

1  1 


12 

I  3 

600 


14 

15 

16 


C2  =  Z  (  M)  *EXPF(  .6405  )  *(  2  795.6*  (  VC  I  (-300.  )-2.  1  If: 5* 

1 V (  I )  )  +  ( 300.- VC  I  (  )  +R2 ( K  )  - V ( I  ( *R  1  ( J ) 

CT2= (-H2  +  S0RTFC ( B2*B2  >-(4.*A2*C2)) )/(?.» A2) 

IFCCT2-I.5)  14,6,6 

6  IFCCT2-2.5)  7,7,14 

7  A3=2.  1  l*Z( M) *EXPF< .8235  1* (  1 .42*R2(K)+530.2*7  CM) *EXPf ( .8235 ) +.25  I  18 
1  *R  1 1  J  )  )  +  1 . 42 *R  1  (  J  )  *R2(  K  )  /  1  .f!3 

E3=2  795.6»Z  C  M ) * EXPF  C . 023  5 ) * ( R  l ( J ) / I . E3  +  2 . 1 1  * Z ( M >  *EXPF ( . 82 35 ) - { 300. 

1  —  V  (  I  )  1/5.2731-RI  (  J  )  *R2  {  K  )  /.I .  L  3 

F3  =  2 . 1  1  *Z  {  M)  *EXPF(  .  82  35  )  Mil  .'42£3*V(  i  )-R2  (K)  )+  1 .42*  (  VC  I  )  *R1  (  J  1-300. 

1  *R2 ( K  )  +V (  I )*R2(K)  ) 

83=F3+F3 

C3=Z  (  M  )  *EXPF  (  .8235  )*(  27*75.6 «  <  VC  I  ) -300.  ) -2.  1  1  F3* 

1  V  (  I  )  )  +  (  300.-VC  I. )  )  *  R  2  (  K  )  —  V  ( >1  )  *  R  I  (J) 

C  T  3=  (  -H3  +  SQRT  F  (  (  B3*B3  )-.{A..'*A3*C3  )  ))/(?.*  A3) 

IFCCT3-I.5)  14,8,8 

8  IFCCTi-2.5)  9,9,  14 

9  A4  =  2.  I1*Z(M)*EXPF(  1 . 00650  * (  1 .42  +  R2IK) +  530. 2  * Z ( M ) *EXPF C I .0065 ) +.25  lid 
1  *R  1  (  J  )  )+1.42»RI(J)*l<2CK)/l.l  3 

84  =  2795. 6»Z(M)*EXPF(  1.00  651  »TrHJ  1/1. E3  +  2.11»Z(M)*EXPF(1.0065)-(  30 
I  O.-VC I  )  )/5.273)-Rl(J)*R2(M/l.fc3 

F4  =  2. 1  l*Z(M)*EXPFC  1.00651  Hi. 42 E 3 «V( 1  1-R2CK)  ) + 1 . 42* ( V (  I 1» 

IR  I  (  J  1-300. *R2( K  )  +V(*I  f +R2C.KH 
B4  =  F  4  +  F4  '' 

C4  =  Z  CM) »EXPF (  1 .0065 )*( 2795. 6*( VI  1 1-300. )-2. 1 1E3* 

1 V  C  I  )  1  +  C  300. -VC  I  )  )  *R  2  C  K  )-VU»  *RJ  (  J) 

CT4=  (  -R4  +  SQRTF  (  (  B4*B4  )  -  (  4  .  * A4  *,C 4  )  )  1/ ( 2. *A4 ) 

IF (C 14- 1.5)  14,  10,  10  ' 

IF1CT4-2.5)  11,  I  1,  1‘4 

CMAA  =  M AX  IFCCT  l,CT2,CT3,CT4 ) 

CMIN=MIN1FCCT I ,CT?,CT3,C 14 ) 

DEL  1 A=CMAA-CMIN 
I F 1  DEL  T  A- .  15)  12,12,  14 

CM AX  1  =  MAX1F(CMAA, CMAX’T  ) 

CMIN1=MIN1FCCMIN,CMINI) 

CONT INUE 

FORMAT! E4.0, 2E 1 5.0, 2T  13.5) 

WRITE  OUTPUT  TAPE  4,600,  V (  1  I  ,  R  I  ( J  )  , R2 ( K  )  , C M AX  1 , CM I N I 
CONTINUE 


CONTINUE 
CONTINUE 
END  FILE 
PAUSE  1 
END 
END 


Fig.  14 
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Fig.  15 


These  wax  limn*  and  •filiUM  neon  currents  do  have  valise  fox  the  analyst 
after  the  design  has  been  chosen,  and  all  that  refine  Is  the  determina¬ 
tion  of  the  particular  maintaining  voltage  in  any  one  circuit*  This  will 
be  of  immediate  use  provided  that  the  neon  tube  voltage  is  close  to  one 
of  the  increments  of  voltage  used  in  the  digital  computer  program;  in 
other  cases,  interpolation  is  required* 

The  digital  computer  program  which  is  flow  charted  in  figure  12  and 
represented  in  FORTRAN  statements  in  figure  13  has  a  computer  running  time 
of  8  minutes  and  30  seconds*  This  output  data  still  requires  interpreta¬ 
tion,  evaluation  and  presentation  in  a  useable  manner  before  it  can  be 
considered  as  a  true  design  asset* 
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Hor  tc  an  evcfcefuti'w  interpretation  of  the  digital,  computer 
printout,  if  ts  i  ««st  advantageous  to  consider  the  nuon  photoconductor 
circuit  or*  sirj  let  has  is*  us  make  the  following  assumptions:  l) 

The  Maintaining  voltage  of  the  neon  tube  Is  fixed  at  60  volts*  2)  The 
knoWi  parameter  has  only  one  value  namely,  2,5,  2}  Temperature  varia¬ 

tions  are  such  chit  the  ambient  temperature  will  vary  between  25°  C  and 
55°  C,  wTith  these  three  assumptions  two  circuits  will  be  considered,  one 
with  compensation  and  one  without  tha  compensating  shunt  photoconductor* 
The  printout  for  the  compensated  and  uncompensated  circuits  are  shown  in 
figure  16,  These  two  printouts  are  of  th*  same  format  as  the  main  program 
except  that  the  maximum  and  uln  imuui  currents  are  now  solely  attributed  to 
temperature  variations* 

A  plot  of  the  design  parameters  R2  versus  Rl,  which  meat  output 
specifications,  namely  peon  current  between  3.5  and  2,5  mi l Hamper as, 
is  constructed  for  both  the  compensated  and  uncompensated  photoconductor 
circuits*  These  plots  are  figures  17  and  18  respectively.  The  uncom¬ 
pensated  photoconductor  circuit  has  n  minimum  value  of  shunt  resistance* 
This  constraint  is  brought  about  by  the  physical  requirement  of  a  suf¬ 
ficient  voltage  drop  across  the  shunt  resistance  to  fire  the  neon  tube. 

The  compensated  photoconductor  circuit  has  a  dark  photoconduct ive  re* 
sistanc*  in  series  with  the  shunt  resistance  prior  to  firing  which  is  in 
chc  rnoguoV.  region.  Therefore,  orfy  the  uncompensated  citcun  requires 
a  miriimcJR  value  of  shuiwt  resistance*  Assuming  a  firing  voltage  of  135 
"  olts  which  is  &n  observed  average  value,  Jnzu  the  following  calcula¬ 


te  ‘Us  for  the  minimum  value  of  R2  at 


2,5  apply s 
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figure  19 
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The  points  i/h^n  plotted  produce  the  dashed  line  in  figure 
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A  }  ■  r  t?  »r  .  €  in  each  p  lot  in  older  t  observe  the 

effect  of  the  change  r  ?•  1  intent  due  only  to  tei  .perature  variations* 

These  numbers  vi thir  area  represent  the  difference  between  the  maximum 

and  minimum  naon  current  in  mil  liamperes  for  that  particular  set  of  design 
parameters  (Rl  and  R2).  A  comparison  of  the  two  plots  reveals  that  compen¬ 
sation  is  actually  achieved  by  insertion  of  the  shunt  photcconductor.  A 
second  conclusion  that  can  readily  be  made  is  that  for  a  particular  value  of 
Rl  in  the  uncompensated  area  plot  the  variation  in  current  is  constant  and 
independent  of  the  value  of  the  shunt  resistance  R2.  This  constant  value 
is  also  the  limiting  value  which  the  compensated  area  plot  approaches  as 
the  shunt  resistance  becomes  infinite.  Analytically  this  last  result 
could  have  been  foreseen  since  for  large  values  cf  shunt  resistance  the 
effect  of  the  compensating  photoconductor  on  the  total  series  resistance 
diminishes  and  eventually  approaches  zero.  This  simple  neon  photoconductor 
circuit  will  be  employed  later  in  this  section  as  an  illustrative  example 
of  best  design  procedure. 

The  digital  computer  printout  of  the  main  program  whether  It  is  the 
hyperbolic  or  "second  order"  approximation  is  of  the  same  format  (See 
Figure  15  for  printout  format).  For  each  value  of  series  resistance  Rl, 
those  values  of  shunt  resistance  F2  which  mark  the  upper  and  lower  limits 
of  successful  circuit  combinations  in  the  printout  are  plotted  on  graph 
paper  producing  areas.  In  the  case  of  tha  hyperbola  only  two  of  these 
plots  are  necessary,  one  for  the  lowest  maintaining  voltage  encountered 

\  T  _ 

and  one  for  the  highest  maintaining  voltage  since  *  -  is  always  nega* 

5  ;  * 

tive.  The  "second  order"  approximation  expression  is  more  involved  and 

it  cannot  be  mathematically  proven  to  have  a  negative  --  £ J  for  all 

( ^  V 

design  sets  to  be  considered.  In  this  case  increments  of  maintaining 
voltage  v;ere  used  ir  the  iigi.tal  computer  program  and  as  many  plots  as 
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volte ^  u  n\,-  '  nt  i*.  ‘  ,  :h  plot  is  ^  $vu *  design  for 

a  particulfr*.  •*  iiAtairln^  iiriu  a  L '  Z*s  and  all  T9s,  In  both 

approximations  the  resulting  areas,  whether  two  at  in  the  case  of  the 
hyperbola  or  many  as  in  the  secorri  order,  can  now  be  overlapped  to  pro¬ 
duce  a  final  area  which  is  common  to  all  of  the  individual  areas.  This 
final  area  if  it  exists  has  satisfied  all  of  the  restrictions  imposed  on 
temperature,  Z,  maintaining  voltage  and  output  neon  current.  The  resultant 
area  for  the  hyperbolic  approximation  is  shown  in  figure  20. 

In  order  to  observe  the  effects  of  more  restrictive  specif ications, 
the  printout  data  of  both  the  hyperbola  and  f second  order*’  approximations 
was  plotted  using  those  upper  and  lower  limits  of  the  shunt  resistance 
that  yielded  an  output  neon  current  between  1.8  and  2.2  mtlliamperes« 

After  having  plotted  the  individual  area  plots*  the  overlapping  procedure 
was  employed  to  find  tha  common  area,  A  common  area  is  non^existan t.  If 
this  specification  of  output  current  had  been  used  vice  1.5  -  2,5  milli- 
amperes,  then  the  absence  of  an  area  indicates  that  there  is  no  design  para* 
meter  combination  that  is  satisfactory  for  the  entire  maintaining  voltage 
range. 

The  amount  of  compensation  effected  can  ha  seen  by  analysis  of  the 
difference  in  maximum  and  minimum  neon  currents  due  only  to  temperature 
variations.  Refer  to  the  simple  example  mentioned  previously  where  the 
maintaining  voltage  was  fixed  at  60  volts  and  parameter  Z  fixed  at  2V5 
(See  figure  17).  The  mapped  points  in  the  interior  of  the  area  plot  in 
figure  17  represent  the  difference  between  the  maximum  and  minimum  neon 
currents  in  mi 1  Hamper es  due  solely  to  temperature  variation  for  that 
particular  set  of  design  parameters  (FI  and  R2),  Inspection  of  the  plot 
reveals  that  the  deviation  due  to  temperature  for  the  same  value  of  R2, 
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b'*  *  i  *  i.. *  ^  largest  r  <  *  .  •  T  « lae  a  I 

stiM  ••  t  #*h«*  ♦  <*  ■  #c  *  ^31  the  s air-e  value  of  *1,  the 

deviation  .*j£  to  ter peiaturfc  h^s  its  silliest  value  vhfc.*?  R2  takes  or.  its 
smallest  j.&ssib]e  valise  inside  the  plotted  area*  The  same  conclusions 
are  true  for  other  similar  plots.  To  see  the  overall  effect  on  the  m  a  in 
program  for  the  hyperbola  and  M second  orderfr  approximations,,  the  difference 
in  neon  currents  due  to  temperature  variations  (called  delta  in  the  computer 
program)  was  restricted  to  first  .015  and  then  .010  milli&ftiperes.  The  re¬ 
sultant  area  plots  for  the  hyperbola  without  delta  and  with  delta  equal  to 
.015  are  shown  in  figure  21.  With  the  program  rerun  using  delta  of  *010, 
a  resultant  area  did  not  exist.  A  comparison  of  the  final  area  curves  with 
and  without  delta  reveals  that  the  left  edge  of  the  resultant  areas  which 
employed  the  delta  restriction  are  moved  further  to  the  right,  where  as 
the  right  edge  is  unaltered. 

The  resultant  area  curves,  since  they  are  in  effect  the  common  area 
of  the  individual  areas*  cannot  be  mapped  for  the  variation  in  neon  currents. 
Any  point  in  the  resultant  area  meets  all  of  the  specifications  that  were 
required  and  is  a  good  design  set.  To  determine  the  best  part  of  the  final 
design  area  we  will  resort  to  inductive  reasoning.  In  inductive  reasoning, 
a  set  of  individual  cases  is  studied  by  the  experimental  method,  and  from 
the  observations  made,  a  general  principle  is  formed.  The  principle  formed 
by  inductive  reasoning  is  reliable  only  when  all  possible  instances  have 
been  examined.  In  the  Individual  areas,  the  mapping  procedure  resulted  in 
the  conclusions  that  the  best  design  set  would  be  the  smallest  value  of  R2 
and  the  largest  value  of  Rl  that  fell  within  the  plotted  area.  Applying  in¬ 
ductive  reasoning  to  this  information,  the  best  set  in  the  resultant  area 
curve  is  also  the  smallest  value  of  P2  and  the  largest  value  of  Rl.  Consider 
at  ion  also  has  to  be  given  to  the  variation  in  the  design  parameters  Rl  and 
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R2 ,  tr  x  i  'uvai^ant  then  the  apex  cf  the  area  would  he  the  best 

design  set.  R1  a*  *?  R2  .  vny  ha/e  tolerances  themselves  of  from  5  to  10  per 
cent  of  their  non  ip  a  1  value,  hence  some  point  that  could  be  the  center  of 
a  circle  whose  radius  is  the  resistance  tolerance  and  where  the  circle  fell 
entirely  within  the  resultant  area  at  the  lowest  part,  would  be  the  best 
design  set. 

The  resultant  area  curves  of  the  hyperbola  and  second  order  approxima- 
tions  are  almost  identical  in  position  and  shape*  Minor  variations  in  shape 
can  be  attributed  to  the  plotting  of  the  individual  computer  incremental  out¬ 
put  data.  As  an  example  of  this  problem  in  plotting  consider  a  design  in¬ 
crement  change  of  2  kohms  in  R1  where  the  value  of  the  design  parameter  re¬ 
sulted  in  a  successful  circuit  prior  to  the  Increment  change,  hut  failed 
after  the  increment  change.  Theoretically  then  the  plotted  curve  could  pass 
through  the  first  point  and  it  also  could  pass  through  any  point  to  the  right 
of  this  original  point  up  to  but  not  including  the  value  of  the  original 
point  plus  the  increment.  Here  then  is  another  factor  in  considering  the 
number  of  design  sets  to  be  investigated  versus  computer  time.  This  factor 
should  not  carry  too  much  weight  since  in  the  final  design,  the  designer 
will  pick  a  point  somewhat  removed  from  the  actual  boundary  curve  itself. 

The  resultant  areas  as  obtained  by  the  two  approximations  should  be  the  same 
for  large  values  of  R2,  The  difference  is  noticeable  where  R2  is  280  kohms 
or  less.  In  this  range  of  R2  compensation  is  more  effective  and  the  second 
order  area  is  a  more  accurate  representation.  However,  in  the  final  analysis 
the  differences  between  the  two  approximations  prove  to  be  very  minor. 

Returning  to  figure  l*5  where  the  voltage  is  fixed  at  60  volts  and  Z  is 
2,5,  it  appears  that  R2  at  R1  s  83  kilohms.  The  digital  computer 

printout  data  when  plotted  reveals  that  the  lower  curve  is  in  effect  the 
locus  of  combinations  of  R]  and  P 2  which  result  in  neon  currents  of  1.5 
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millicv]  5  •  i .  .i  il  ar  locus  uOi  cuiterts  of  2*  > 

millan.peres*  These  can  j?  w  ctuaV.y  only  approximate  locii  of  these 

constant  current  values  since  the  joints  which  were  used  to  construct  the 

curves  were  obtained  by  using  incxorc&ntal  values  of  design  parameters,, 

2 

The  system  equation  which  was  of  the  form  AI  t  BT  +  C  »  0  was 
solved  explicitly  on  the  digital  computer  for  as  a  function  of  five 
parameteres „  If  the  system  equation  is  rewritten  so  as  to  express  R2  as 
a  function  of  the  remaining  parameteres,  the  result  using  the  hyperbola 
approximation  would  be: 


Using  a  neon  current  of  2*5  mil  liamperes,  a  voltage  of  60  volts,  a 
Z  of  2.5  and  a  mean  temperature  of  40 °C,  the  value  of  El  to  make  the  de¬ 
nominator  of  the  equation  zero  is  85  kohms .  Using  this  equation  similar 
results  could  be  predicted  for  other  maintaining  voltages,  namely  that  in 
the  plot  of  design  parameter  R2  versus  design  parameter  Hi  there  are  values 
of  Hi  where  R2  goes  to  infinity.  Since  the  resultant  area  is  constructed 
from  the  individual  area  plots  we  cannot  expect  to  have  a  closed  resultant 
design  area  plot. 

The  insight  gained  from  the  interpretation  of  the  digital  computer  data 
for  the  circuit  synthesis  problem  revealed  that  the  constructed  boundary  line 
are  in  essence  locii  of  constant  neon  currents,  namely  the  extremes  of  the 
specified  neon  current  limits  of  1.5  and  2,5  milliamperes.  This  suggests  an 
alternate  technique  for  solving  the  circuit  synthesis  problem.  This  alter- 
nate  technique  lends  itself  to  solution  using  an  analogue  computer. 
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7.  De s.fijv.  ^ -it  ionic  \Pi 1  ogue  Computer  Technique 

The  technique  t o  be  employed  with  the  analogue  corputer  parallels  the 
development  of  the  digital  computer  technique  to  a  certain  extent.  Known 
and  design  parameters  are  as  defined  previously.  The  resultant  circuit 
equation  is  manipulated  to  express  cne  design  parameter  explicitly  as  a 
function  of  the  remaining  design  and  knovm  parameters*  Increments  of  the 
known  parameters  and  end  points  of  the  known  parameters  a  re  determined  as 
in  the  digital  computer  technique.  The  explicit  circuit  equation  is  then 
converted  into  an  electronic  analogue.  The  electronic  analogue  computer  is 
set  up  so  that  its  equations  have  the  same  mathematical  form  as  the  equation 
for  the  physical  system  under  considerat ion.  In  the  computer,  voltages  are 
used  to  represent  various  physical  quantities  such  as  velocity,  position,  etc. 
Scaling  factors  are  used  to  convert  the  physical  quantities  into  appropriate 
numerical  values  in  the  computer.  The  mathematical  functions  of  addition, 
multiplication,  etc.,  are  performed  by  electronic  operational  amplifiers 
using  scaled  voltages  as  inputs. 

large  variations  in  knovTi  parameters  can  be  readily  undertaken  provid¬ 
ed  the  range  in  variation  is  considered  during  the  scaling  process.  Incre¬ 
mental  charges  in  the  knoxro  parameters  are  accomplished  by  merely  varying 
the  coefficient:  potentiometers  that  are  associated  with  the  parameter  under 
invest igation.  Radical  changes  in  the  known  parameters  may  have  to  be 
handled  by  varying  the  components  of  the  operational  amplifiers,  or  by  the 
insertion  of  one  or  more  additonal  amplifiers. 

A  design  parameter  that  is  implicit  in  the  circuit  equation  is  varied 
by  the  application  of  a  ramp  voltage  input  causing  the  explicit  parameter 
output  voltage  to  be  the  continuous  solution  for  a  fixed  set  of  known  para¬ 
meters.  Utilization  of  an  X-Y  plotter  will  present  a  graphical  representation 
of  a  design  parameter  voltage  as  a  function  of  another  design  parameter  voltage. 
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Eich  tii  o  that  t.  !i  .<  <ii  ;  ■  •  •  ,  is  *oj  <<  t  curve  ul!l  ’  -»  tra-.  cl  «.  at  on 

tHs  X"  x  plotter ,  -hu  ut  b s «'  u  ' . r. c~  n  parameters  that  cor*  he  varied  .oti  any 
one  sheet  of  the  i-i  p  .otter  hecooes  oviir.et  v»liifi  .5  loss  of  identity  in  the 
plotted  curves  results.  This  pro&dure  is  repented  until  graphical  plots 
are  obtained  ter  every  increment  of  each  of  the  known  parameters.  These  plots 
are  antiscaled  on  the  X-Y  plotter  so  that  the  graphic®!  plots  which  are  ob¬ 
tained  in  this  manner  are  design  parameters  as  used  in  the  original  physical 
system  prior  to  the  analogue  conversion.  Each  plot  has  only  investigated  one 
or  more  variations  in  known  parameters,  tf  all  of  the  plots  are  overlayed, 
the  area  which  is  common  to  the  acceptable  area  of  ell  the  graphs  is  the  final 
design  area  desired.  Any  point  in  this  final  common  area  meets  all  of  the  re¬ 
quired  specifications. 

Rearrangement  of  the  basic  circuit  equation  employing  the  hyperbola 
approximation  for  the  photoconductors  resulted  in  the  form: 


2* 
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where  B  3.W  £ 

The  neon  current  I  is  now  an  end  point  type  of  known  parameter  rather 
than  the  output  quantity  as  in  the  digital  computer  technique.  The  two  neon 


current  values  of  interest  are  1.5  and  ?,5  milliairperes.  Physical  quanti¬ 
ties  such  as  ohms,  amperes,  and  maintaining  voltage  in  the  circuit  equation 
are  transformed  intc  voltages  for  use  in  the  electronic  analogue  by  appropri¬ 
ate  scaling.  Known  parameters  of  temperature  and  Z  are  divided  into  in¬ 
crements  as  before.  A  change  in  the  value  of  any  known  parameter  is  accomp¬ 
lished  by  altering  the  individual  potentiometers  that  are  associated  with 
the  parameter  under  investigation.  Since  the  range  of  all  the  known  para¬ 
meters  vary  by  less  than  ?  factor  of  ten,  it  is  possible  to  utilize  one  analo¬ 
gue  computer  circuit  and  account  for  all  of  the  variations  by  suitable  scal¬ 
ing  of  the  physical  parameters.  .  schematic  of  the  analogue  computer 
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Mult. 

2BZI(R1  +  103BZ) 

-2BZ(300-V) 


amplifier  **  j  7 1  .-st t ated  u  figure  22, 

A  step  voltage  a ppliei  to  iniergrat ing  amplifier  F  produces  a 
ramp  voltage  which  corresponds  to  design  parameter  Rl.  This  continuous  vari- 
ation  of  design  parameter  Rl,  rith  a  fixed  set  of  temperature,  Z,  maintain¬ 
ing  voltage  and  neon  current  will  produce  a  continuous  solution  for  the  ex- 
plicit  design  parameter  R2.  Potentiometers  a4,  a5,  a7,  alC,  nil  and  al2  are 
held  constant  at  one  value  throughout  the  problem.  These  constant  potentio¬ 
meter  settings  are  determined  by  scaling  physical  quantities  which  are  in¬ 
variant,  The  setting  of  the  variable  potentiometers  at  any  one  instance 
depend  upon  the  known  parameter  that  is  being  investigated.  Table  #2  lists 
the  dependence  of  each  known  parameter  on  its  associated  potentiometers. 


Table  ^2 


Variation 

in 

I 

Potentiometers 

a3. 

a8s  a!5,  nl6 

Variation 

in 

T 

Potentiometer 

a2 

Variation 

in 

2 

Potentiometer 

al 

Variation 

in 

V 

Potentiometers 

a9. 

al3,  a 14 

Figure  23 

Each  of  the  known  parameters  one  at  a  time  is  systematically  varied  by 
changing  their  associated  potentiometers.  The  tv;o  voltages  in  the  electron¬ 
ic  analogue  which  represent  design  parameters  Rl  and  R2  are  automatically 
plotted  by  inserting  these  continually  varying  voltages  as  inputs  to  an 
X-Y  plotter  during  any  one  solution  which  corresponds  to  a  particular  knov/n 
set  of  parameters.  The  output  of  the  X-Y  plotter  is  a  graphical  representa¬ 
tion  of  design  parameter  R2  versus  design  parameter  Rl.  A  typical  output 
curve  for  Z  (photoconductor  resistance  magnitude  multiplying  factor)  corres¬ 
ponding  to  .4  appears  as  in  figure  24,  This  curve  includes  the  effect  of 
temperature,  maintaining  voltage  and  neon  currents.  Solid  lines  represent 
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maintaining  voltages  of  volts  rt  \  dashed  lines  represent  maintaining 
voltages  of  ICO  volts.  The  three  temperatures  investigated  are  25°C,  4Q°C, 
and  55 °C  which  are  recorded  from  right  to  left  in  groups  of  three  lines. 

For  each  maintaining  voltage*  the  lover  gr< up  corresponds  to  1.5  mil liamperes 
and  the  remaining  group  corresponds  to  2,5  milliamperes.  Similar  curves  for 
other  values  of  Z  appear  in  figures  25,  26  and  27, 

As  was  mentioned  previously  in  section  four,  the  known  parameters, 
temperature  (T),  and  photoconductor  resistance  magnitude  multiplying  fac¬ 
tor  (Z),  were  established  as  incremental  types  of  parameters  since  ---■  - — 


and 


yr 

j  * 


could  not  be  proven  to  be  of  constant  sign  for  the  entire  range 


of  variation,  Examination  of  figures  24  through  27  discloses  that  for  the 
maintaining  voltage  equal  to  60  volts,  the  three  temperature  curves  of  the 
lower  group  (neon  current  equal  to  1,5  mliiamperes)  cross  over  each  other. 


This  cross  over  is  indicative  of  a  point  of  inflection  in  or  stated 

J  1 

in  other  terns  the  value  of  L  changed  sign  at  the  point  of  temperature 

J  T 

cross  over,  By  overlaying  figures  24  and  25  a  similar  effect  would  be  noticed 

for  Z,  namely  that  for  one  value  of  temperature,  maintaining  voltage  and  neon 

current  the  overlayed  curves  corresponding  to  two  different  values  of  Z 

would  cross  over.  This  cross  cvei  indicates  where  a  point  of  inflection  ex- 
\  T 

ists  for  —  ,  An  overlapping  technique  as  applied  to  the  digital  computer 

)  r* 

<y 

areas  is  used  in  finding  the  area  common  to  the  individual  areas,  re¬ 

sultant  design  area  is  figure  28.  Since  all  of  the  points  of  inflection  for 
temperature  and  photoconductor  resistance  magnitude  multiplying  factor  lie 
outside  the  resultant  design  area,  these  quantities  are  in  effect  end  point 
type  of  known  parameters  if  the  term  as  used  here  applies  only  to  the  result¬ 
ant  design  area. 

The  design  area  as  obtained  from  the  analogue  computer  in  conjunction 
with  the  X-7  plotter  meets  the  requirements  of  circuit  performance  with  any 
set  of  known  parameters.  Thi  Hues  forming  the  boundaries  of  the  area  are 


",  ,  j  -  n-  in.  M  •  .  i  (,  S' 


O  V 


/3o 


ejw,cc  l  s  *  iv>  .  bat  •  i«  t  i  - 

ital  ■  pit  i*st  :e  t  t^-eie  r  combination  -u*r  -t 

ed  from  this  design  area  as  was  obtained  from  the  digital  cornpj 
there  is  no  add  it.  nal  information  available  about  the  resultant 
upon  which  the  designer  can  base  a  choice*  This  perhaps  Is  its 
as  compared  to  the  digital  computer  technique* 


be  determin¬ 
er  d/'ta*  since 
design  area 
sole  weakness 
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R2  IN  //(LOH/nS 


C  *  1  as  ...  'ns  * 

The  profc len  cf  how  lo  a  circuit  to  meet  specifications  h  le 

using  r aremefct- vs  which  have  w.fde  ranges  of  variation  can  be  solved  using 
either  of  the  computer  techniques  described*  The  designer  has  merely  to 
select  design  parameters  from  a  final  area  with  the  assurance  that  the 
circuit  will  meet  the  output  specifications* 

Both  the  digital  and  the  analogue  computer  techniques  provide  adequate 
design  data*  Each  method  has  advantages:  there  is  little  difference  in  the 
amount  of  time  involved  from  start  of  problem  until  solution*  This  time 
naturally  depends  upon  the  designer’s  experience  and  familiarity  with  a 
particular  computer*  The  accuracy  is  comparable  in  either  method;  while 
the  digital  computer  is  more  accurate  for  any  individual  calculation*  the 
analogue  has  the  advantage  of  providing  a  continuous  solution  rather  than 
an  incremental  one.  The  analogue  computer,  when  used  with  an  X-Y  plotter 
provides  the  area  curves  directly  without  the  plotting  necessary  to  reduce 
the  digital  data*  However,  the  digital  computer  can  be  easily  programmed 
to  provide  addilion&l  data  useful  in  analysis  of  a  specific,  designed  cir¬ 
cuit*  There  were  no  logic  circuits  employed  in  the  analogue  solution, 
therefore  the  significant  test  for  delta  I  was  not  made  in  this  solution* 

For  the  specific  circuit  used  to  demonstrate  these  techniques,  it 
should  be  noted  that  the  scheme  for  temperature  compensation  did  achieve 
the  desired  compensation  as  reference  to  figures  17  and  18  will  show.  The 
additional  labor  involved  in  employing  the  more  complex  second  order  approxi¬ 
mation  is  not  warranted  since  there  is  no  significant  difference  between  the 
final  design  areas  obtained  with  the  two  approximations*  Figure  29  compares 
the  results  obtained  with  the  two  approximations  using  the  digital  computer 
and  the  hyperbolic  approximat ion  using  the  analogue  computer. 
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Additional  investigation  should  be  conducted  to  extend  the  study  to 
include  an  analysis  of  the  circuit* s  transient  behavior*  A  method  of  in¬ 
cluding  logic  circuitry  in  the  analogue  circuit  would  permit  the  additional 
test  for  delta  I  to  be  conducted*  A  circuit  involving  the  selection  of 
three  design  parameters  presents  an  interesting  problem*  In  this  case, 
instead  of  acceptable  area  graphs,  the  result  would  be  a  three  dimension¬ 
al  acceptable  volume  plot  of  the  design  parameters. 
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APPENDIX  I 

MAT  (EMATICAL  DERIVATION  LEADING  TO  SECOND  ORDER  EXPRESSION  FOP 
P  TOTOCONDUCTOR  RESISTANCE-NEON  CURRENT  RELATIONSHIP 

To  obtain  a  mathematical  expression  which  approximates  the  photoconductor 
resistance  variation  with  neon  current  shown  in  curve  B,  the  relationship 
is  assumed  to  be  of  the  form: 

J\  P,,L  t  B  1  t  ('  K'j,  I  i- 

The  following  points,  selected  from  the  curve  were  used  to  determine  the 
unknown  coefficients  A,  B,  C  and  D. 


I  in  ma.  1,50  1.80  2,20  2.50 

R  in  Kohm  6.71  5.10  3.88  3.37 


A  (  L  n  i(j  5  )  t  I  i  (  '  2  XJ  o  ')  f  (  (  i  "J  •  p 
A  ( i-  0  (j  I'  )  1  B  (  I  i  /  A:  '  )  t  (  o  ;  //<;"  J  *  ' 


(1) 


'  (2) 


/\  (j  ?  -V  5  )  t  B  (jl.  I  \/c  ])  i  C  (  h  4**  X'  ■  ? 


(3) 


/\  ^  )f  ?  ;  f  )  f  h  >  X  .  j  r  v  <_  ’/A 


(4) 


ft'*’  A  -  *  <  •  ''  B  c  /.  LI  Af  t  5  C  -  C 


(1)  -  (2) 


j  |  I  A  3  i<  /  C.  J3  f  5  /  x  /  • 


(3)  -  (4) 
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Subtracting' 


/  -S 


1  .. 

f-  <  <  A 


•V  < 


(3a) 


i  I  ,  ci  i  x 


v  *  »• 


•f  y  •  x  i 


2,  7  X ;  a  3 


2  x  !  c  /3 


Ji  ■< 


y  ,7?  y  s 


<  / 


(4®) 


Subtracting, 

2xio  rS  x 

fi  ll  xjo" 

Substituting  In  (4a)  and  simplifying. 


X,  6  '  (  M 

u.j*>  X  / 


I  V  .■  ./*  ^  i* 


J  *>* 


J-  **  ;  X 


P  i  A  X 


for  t''  / 

/^i/  i  X  /  <-  *  ;  f  j 


C  L 


T 


J_  .<  ‘  i  -A  /  : 

h  'it  1  X  /  iJ  jT  I 


A  J  -X  / 


i  *  2. _ t 

l  H  •  ‘1  **.s:  f 


To  incorporate  the  effect  of  temperature,  rewrite  in  the  following  forms 


A 


;  n  .L 


A 


V  i  *  1  ) 


Jf 
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*  final  f c rni  of  e  sjno.  '  *u  in l*  >t  toxi  atio  _' s 


l; 


■  ■  ! 


!  ■;  -  ■  -  ; 
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Heat  ion. 
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